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Abstract 

We discuss the conditions for which the non-equilibrium work relation is valid by 
means of thermodynamic and microscopic arguments. 
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1 Introduction 



The question of dissipation and the approach to equilibrium in mesoscopic and 
macroscopic systems has been a long standing problem since the early years 
of statistical mechanics as a science. Already at the beginning of the XX cen- 
tury Einstein and Ehrenfest were unsatisfied with the fundamental principles 
of statistical mechanics postulated by Bolztmann and Gibbs, since in their 
opinion those principles lack of a sound microdynamical basis. (1; 2; 3) Their 
main idea was that the real problem concerns the approach towards the equi- 
librium state, and not equilibrium in itself. So, the knowledge of the dynamics 
of the irreversibility and the approach to equilibrium that follows enables one 
to draw the framework embodying equilibrium as a limiting behavior. 
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This subject has acquired a renewed interest, mainly due to the technological 
advances in biophysics and microrheology that make a more detailed analysis 
of the question possible. (4; 5) 

Related to this problem, attempts have been made to go beyond equilibrium 
in order to establish some rules for extracting information from irreversible 
processes. (6; 7; 8; 9; 10; 11) Recently it has been claimed that there are only 
a few known relations in statistical mechanics that are valid for systems arbi- 
trarily driven far from equilibrium. One of these relations is the one proposed 
by Jarzynski (12; 11) for which equilibrium properties can be obtained from 
non-equilibrium work measurements. For Hamiltonian isolated systems and 
systems in contact with a heat bath through weak interactions, theoretical 
justification for this assertion has been given in terms of ensembles of tra- 
jectories represented by means of a phase space density. (11) However, from 
the thermodynamic point of view this assertion seems to be true only for 
very specific situations. (13) As a consequence, it is important to stipulate the 
conditions in which one can make use of this relation. On the other hand, 
the knowledge of the range of applicability of this and other non-equilibrium 
relations is deeply rooted in the understanding of the microscopic bases of 
irreversibility. (6; 7; 8; 9; 10) 

In this article, we will use thermodynamic and statistical arguments to show 
that the non-equilibrium work relation is valid only in isothermal and near- 
equilibrium conditions. Moreover, we will revise the validity of the original 
derivation of this non-equilibrium work relation (11) and then establish a 
quantitative criterion for determining these conditions by making use of a mi- 
croscopic theory involving the BBGKY (Bogoliubov-Born-Green-Kirkwood- 
Yvon) hierarchy of equations (14). 

The paper is organized as follows. In Section 2, we analyze the non-equilibrium 
work relation by considering the thermodynamic definition of the minimum 
work done on a system to change its state. After this, in Sec. 3 we analyze the 
problem of irreversibility in the context of a microscopic theory and propose 
the entropic time as a quantitative criterion to distinguish quasistatic and 
nonquasistatic processes. Finally, the last section is devoted to summarize 
and discuss our main results. 



2 Thermodynamics and cumulants 

In the very general system driven out of equilibrium does not satisfy 

the condition of thermal equilibrium with the bath, since the temperature T of 
the system differs from that of the bath T (15). Among others, typical systems 
illustrating this situation are glasses and supercooled colloidal fluids. (16; 17) 
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Since irreversible processes are present in this general case, the total work R 
performed on the system in order to change its state satisfies the relation (18) 

R > AE - TnAS, (I) 



where AE and AS* are the energy and the entropy variation of the system in 
the process. Hence, in the case of a reversible process, from Eq. (1) one can 
infer that the minimum amount of work R m in necessary to reverse the state 
of the system is given by 

R min = AE-T AS. (2) 



In the isothermal case, when the system is in contact with a heat bath, T = T , 
Eq. (2) can be used to define the variation of the Helmholtz free energy of the 
system AE 

R mm = A(E - T S) = AE (3) 



Therefore, thermodynamics requires that the change of the free energy of the 
system be equal to the work performed on it only in a reversible process 
in which at least the initial and final states satisfy the thermal equilibrium 
condition T = T . 

In Ref. (11), the following equation has been established 

AE = -k B T \n(e- R/kBTo ), (4) 



where ks is the Boltzmann constant and () represents the average over an 
ensemble of measurements of the work R. 

In view of the thermodynamic relations (l)-(3), one may legitimately ask 
about the conditions of validity of Eq. (4). To answer this question, it is 
very important to take into account that, in general, an arbitrary external 
perturbation will produce an internal irreversible process in which the system 
approaches equilibrium, in accordance with the Le Chatelier principle. (18) 

For simplicity we will assume that work is done in isothermal conditions; 
therefore, we have 

R = AF + R diss , (5) 



where Rdiss is the dissipated work. (18) Now, by taking the exponential of Eq. 
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(5) one simply obtains 

e R/k B T _ e AF/k B To e R Mss /k B T 



(6) 



After averaging this expression over the ensemble of measurements of R one 
obtains 

/ e -R/k B T \ _ e -AF/k B T / e -R dias /k B T \ (j) 



where we have taken into account that AF is a constant quantity between 
two equilibrium states. Therefore, from Eq. (7) it follows that the relation 
^ e -R/k B T \^ ^ e -AF/k B T j g va ^(j on iy jf the time scales characterizing the vari- 
ation of the external parameters are larger than the time scales characterizing 
the decay of the irreversible processes taking place within the system, since 
then Rdiss ^ AF. A quantitative criterion for this time scale will be given in 
the following section, where these irreversible processes will be analyzed on 
microscopic basis. 

For situations not far from equilibrium the implications of Eq. (6) can be 
viewed in the context of the linear response theory. Thus, by expanding the 
exponential containing Rdi SS and averaging over the ensemble of measurements 
of R. Up to first order, (15) the result is 

R/k B T \ ^ AF/k B T 



{Rdiss) 

k B T 



0(2) 



which seems to imply that Eq. (4) does not contain the first-order correction 
related to the response function of the system, included here in Rdiss- 

This conclusion is in accordance with the requirement that Eq. (4) be valid 
for fluctuations of the work smaller than ksT as assumed in Ref.(ll): "the 
fluctuations in the work R must not be much greater than /c^T, if we are to 
have any hope of verifying Eq. (4) experimentally...". A similar result can be 
obtained, for example, by expanding the average (e R / fcsTo ) around the mean 
value (R) assuming that the distribution of work is a Gaussian. Due to the 
symmetry of this distribution one obtains 

i p R/k B T \ _ JR)/k B T 



1 + 



(R 2 ) - W 

2(fc B T )2 



(9) 



which as compared with Eq. (8), leads to the estimate 

(Rdiss) {R 2 ) — (R) 2 

IbYo 2(k B T )i ' 



(10) 
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This is precisely the result obtained in Ref. (11) in the paragraph that follows 
Eq. (12) in this reference. In addition, given that (R) = AF + (Rdiss), by sub- 
tracting (8) from (9) and taking into account (10), one obtains e ( ~ Rd,ss ^ kBT ° = 
1, which means that (Rdiss) Ab^o ~ 0. Our linear response analysis, although 
it may seem naive, runs parallel to the one underlying in the assumptions con- 
cerning the work fluctuations in Ref. (11). A more elaborated linear response 
analysis has been done in Ref. (19), in which the authors conclude that what 
the non- equilibrium work relation really offers is not the thermodynamic free 
energy difference but merely an upper bound: "most significantly, AAj is not 

the free energy change predicted by thermodynamics ". This assertion is 

relevant since it also indicates that Eq. (4) gives only an approximate value of 
AF, thus implying that this is not valid in general, contrary to what is com- 
monly claimed. This is a widely extended misinterpretation which we attempt 
to clarify here. 

The previous results highlight the thermodynamic implications of the two 
assumptions underlying Eq. (4). However, a deeper analysis can be performed 
by noticing that l^e R / kBT °) is the moment-generating function of the probability 
distribution of R. (20; 21) Hence, in a general case we get more clarity if we 
perform a series expansion of the quantity — fc^To ln(e~ i? / fcsT °), 

t^i n] - (-Wo) 
where K n is the n-th order cumulant. 

Using now the theory of Thiele semi-invariants, (20; 22) we can define the n-th 
order cumulant in terms of the moments of (-R fc )- F° r the first few values of 
n, we have 



ki = {R) (12) 
K 2 = (R 2 )-(R) 2 (13) 

k 3 = (r 3 )-3(R 2 )(R) + 2(R) (14) 

K4 = _ 4 (i? 3 ) (R) - 3 (i? 2 ) + 12(i? 2 ) (R) — 6 (R). (15) 

It can be shown(21) that for a Gaussian distribution, K n = for n > 2. So, 
the value of the quantity — ksT h\(e~ R ' kBT °) depends on the probability dis- 
tribution, and therefore cannot coincide with the thermodynamic free energy 
since the thermodynamic free energy must be independent from the statis- 
tics. As we have shown above, this coincidence occurs only in the Gaussian 
case when the system is not arbitrarily far away from equilibrium {i.e., in the 
fluctuation-dissipation regime). A cumulant analysis has been also performed 
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in several papers, in Ref. (11) itself or in Refs. (19; 23) just to cite a few 
exemples. Nonetheless, the approach in these works differs from ours since in 
these the non-equilibrium work relation is taken for granted and the cumu- 
lant expansion serves to test how great the statistic must be to compute the 
free energy. For us however, the cumulant expansion is less restrictive since 
the only thing that this gives us for sure is information about the probability 
distribution of work fluctuations and very little else. 

Let us illustrate our point with a simple example which shows that Eq. (4) 
is not general. A typical distribution in non-equilibrium situations is the chi- 
squared or gamma distribution, (24) for which 



- k B T ln(e- R / kBTo ) = -k B T \n(e-W AF ^ AF / kBT ^) = (16) 

where work has been measured in units of AF. The sum in Eq. (17) is a 
Mercator series which is convergent when 2AF/k B T < 1. Hence in this case, 
we obtain 



k B T Q In(e-^-o) = -M" M In ( 1 - 2^- ) * (11 



2 AF V k B T 
(R) = (Rdiss) + AF = AF (\ + + AF (19) 

Therefore, clearly the non-equilibrium work relation is not satisfied in this 
case. Thus, by considering nonequilibrum situations where the distribution of 
Rdiss can be fitted to a chi-square distribution 

F(x)dx = l/[2r(//2)] (x/2) f/2 ~ 1 e~ x / 2 dx, (20) 



with variance a 2 = 2f and mean m x = f, with x = Rdi SS /k B T Q and / ~ (x), 
we can compute 

oo 

= J e- x F(x)dx = (l/3) //2 ± 1 (21) 



Hencen, since according to Eq. (21) (e~ x ) ^ 1 one concludes that (e~ R l kBTa ) ^ 
e -AF/k B T w \yi c \y violates the conditions we have found through Eqs. (8)-(10). 
These results clearly show that unlike what is claimed in the current literature 
(12), for arbitrary large fluctuations the non-equilibrium work relation is not 
applicable and thus is only valid for equilibrium fluctuations. 
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Experimentally, the mean value ( R) is obtained through an arithmetic mean. 
A Gaussian distribution centered in (R) is used to account for data dispersion. 
Since the number of experiments performed is finite, this distribution corre- 
sponds to the statistical distribution of frequencies of the measured values 
of R. However, as is emphasized in Fig. 1, this statistical distribution does 
not necessarily coincide with the limiting distribution whose maximum value 
corresponds to the free energy difference AF. As a consequence, it is quite 
possible that the measurements will have values near those of the limiting 
distribution, and any difference will be interpreted as a statistical error and 
not due to other kind of bias. 

To conclude our analysis of Eq. (4), we will consider a system far from equi- 
librium which is in contact with a heat bath, whose entirety constitutes an 
isolated system(ll). Under these conditions, it is very important to note that 
the temperature T oc P~ x is a function of the energy, (25; 26) regardless of 
the size of the system (27). There is also more recent literature emphasizing 
this fact, (28; 29; 30). Therefore, for a time-dependent Hamiltonian system, 
T is not a thermodynamic variable in itself but rather a time-dependent pa- 
rameter. Hence, the use of a constant equilibrium temperature To oc Pq 1 as 
a prefactor in the exponentials of Eqs. (7) and (8) of Ref. (11) is not appro- 
priate, since its time dependence must be taken into account. In our opinion 
this non-adequate treatment of the temperature is also present in Ref. (31). It 
should be noted that our perception about the temperature is coincident with 
the same objection raised in Refs. (13; 32). In view of this, we may follow Ref. 
(11) in order to write 



/(*; t) e -* u '( 2; *) = £-i e -/W*D) e -A>«(*;t) (22) 
= Z - 1 e- /3o ^ (z;t) , (23) 

where w(z; t) = H\(z; t) — Hq{zq) is the work done along a trajectory z(zq, t) 
with initial condition Zq. However, since for a finite switching time t s (3{t s ) = 
Pi 7^ fa, one obtains 

J dzexp[-p Q H x ] ±Z X = J dzexpl-PiH^, (24) 



in contradiction with the result given in Ref. (11). Additionally, it must be em- 
phasized that the definition of work w in Eq. (23) is not consistent with the one 
corresponding to thermodynamic processes. The consistency between these 
two definitions arises only when the processes considered are slow enough, 
i.e., for adiabatic processes for which the entropy is a constant and therefore 
constitute reversible processes. (18) 

As a consequence of the previous analysis, it follows that it is not correct to 
identify the quantity w(ll) with the thermodynamic work R which, through 
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Eqs. (3) and (5), constitutes the definition of the free energy difference AF. 
Thus, after using Eq. (7), one can write 



" ( e -R/k B T \ - 



which is the general relation arising from thermodynamic arguments. 



3 Microscopic analysis: The entropic time 



The discussion in the previous section allows us to conclude that, far away from 
equilibrium, the approximations made in order to derive non-equilibrium work 
relation's (4) are no longer valid. Hence, the evaluation of the dissipated work 
becomes necessary. 

In order to achieve this objective and to obtain a quantitative criterion estab- 
lishing when Eq. (4) can be applied, we will write the work dissipated by a 
system in contact with a heat bath and due to the action of an external force 
Xj in the form(33) 

Rdiss = XjA^j, (26) 

where £j is an internal parameter characterizing the state of the system. Since 
in view of Eq. (26) Rdiss is related to the entropy production of the system 
<t, we may write Eq. (7) 



l e -R/k B T ^ ^ e -AF/k B T 



^ + 0(2) 



(27) 



where r s is an elapsed time. From this equation it follows that it is necessary to 
estimate the characteristic relaxation time of the entropy production in order 
to determine when Eq. (4) is valid. 

This task can be accomplished by means of a microscopic theory that we 
will summarize here. For an isolated N-body system, in Ref. (34; 35; 36) it 
was postulated the non-equilibrium entropy S in the BBGKY scenario as a 
functional of the set of s-particle reduced distribution functions (s < N), 
represented in the distribution vector f 



S = -k B tr {f In (f^f ) } + S eq = (2* 
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N 1 f fs 

~ k B^2 — / f s ln—^-dx 1 ...dx s + S eq , (29) 
s=i s - J Jeq,s 

which generalizes the Gibbs entropy postulate(7; 8; 9). In this relation S eq 
is the (thermodynamic) entropy at equilibrium and f eq is the distribution 
vector that corresponds to the equilibrium state. Hence, f eq is the equilibrium 
solution of the BBGKY hierarchy, in other words this is a solution of the YBG 
(Yvon-Born- Green) hierarchy (14; 37), and therefore is not a Boltzmann-like 
function. Additionally, the non-equilibrium entropy (28) reaches its maximum 
value at equilibrium, when f = f eq . 

In the framework of the BBGKY description, the system is assumed to be 
a mixture of s-particle interacting fluids in the phase space and each fluid is 
made up of particle-clusters of equal size. (34) As a consequence of the inter- 
action among those fluids a compressible multiphase flow is established in the 
phase space, so it has been proved (34) that Eq. (28) is a monotonically in- 
creasing function in time that properly describes the regression to equilibrium 
of a system originally under non-equilibrium conditions. 

In Ref. (34) , it was also shown that the dynamics of the non-equilibrium 
entropy (28) follows from the dynamics of the distribution vector f given 
through 

^ = £f(t), (30) 



which is the generalized Liouville equation expressing the BBGKY hierarchy 
in a compact way, where £ is the generalized Liouvillian. As it is known, Eq. 
(30) is obtained by projecting adequately the Liouville equation onto each 
one of the s-particle phase space. Therefore, using Eqs. (28) and (30) one can 
compute the entropy production 




where AjF^ is the averaged non-equilibrium force between the j-th particle 
and the particles of the remaining fluids. Moreover, Jj arises from deviations 
of the distribution function with respect to equilibrium. 

From Eq. (31) we may estimate the relaxation time of the irreversible processes 
taking place in the system. To this end, note that the product Jjdx\...dx s 
has dimensions of velocity whereas AjF has dimensions of force. Introducing 
the characteristic velocity v and taking into account that AjF contains the 
interaction forces among the components of the system, then AjF ~ 0o/ r o 
with 0o an d ?"o being the characteristic energy and the characteristic length 
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of the interaction potential. Thus, the entropic time scale t s = mvro/4>o (with 
m the mass of the system) related to the approach to equilibrium is the time 
associated to the relaxation of the non-equilibrium entropy to its equilibrium 



value. If we assume that v ~ y^sT/m with m the mass of the system, we 
finally obtain 



which has the expected physical limits at low and high temperatures. Eq. 
(32) implies that experimental measurements can be interpreted in the con- 
text of Eq. (4) only if the time elapsed before taking a measurement in each 
point of the trajectory is larger than t s . It is important to note that the de- 
pendence on mass of Eq. (32) implies that the time of decay of irreversible 
processes in nanoscale systems could be very short, thus making them ap- 
propriate candidates to make the misleading interpretation that Eq. (4) is a 
far-from equilibrium relation. 

Here, it is convenient to consider that the thermodynamics of small systems, (27) 
introduces corrections to the thermodynamic variables and state functions 
which make them different from their corresponding macroscopic counterparts. 
The magnitude of these corrections depend directly on the size of the ensem- 
ble of systems which is considered. (27) From a dynamical point of view, for 
small systems these deviations from the macroscopic behavior are not small 
and may affect, in general, not only the extensive quantities but also the in- 
tensive ones (as, for example, the Massieu function). This fact is important, 
since these fluctuations modify the dynamics in such a way that may introduce 
multiplicative stochastic noises in the corresponding evolution equations for 
the variables determining the state of the system. Often, these multiplicative 
fluctuations lead to polynomial corrections to the local Gaussian distributions 
associated with the fluctuating quantities. As shown in the previous section, 
these corrections will affect the cumulant expansion of the corresponding dis- 
tribution. The assumption that for arbitrary nonequilibrium processes the 
higher order terms of Eqs. (8)-(13) cancel between them in general, seems to 
be incorrect. 

In the previous paragraph we analyzed the irreversible behavior of an isolated 
system. However, since most of the experiments are performed in the presence 
of a bath, we should generalize our analysis. 

In the case when the system is in contact with a bath, it is necessary to 
extend the full phase-space in order to account for the presence of the bath. 
In this framework, the volume elements of the new phase space are conserved 
and therefore the Liouville theorem is satisfied. Thus, we will include a set 
of pairs of conjugated generalized coordinates (A, A) so that the Hamiltonian 
is now a function of A and A, H (p, q, A, A). (18) Here, A represents the set of 





(32) 
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external parameters determining the state of the system and accounting for 
the interaction with the bath and A the set of conjugated generalized momenta 
(which can be also understood as external currents). 

Therefore, in this case equation (30) is rewritten in the form 



where d\j dt is the convective derivative associated with the external currents. 
Following the procedure that bring us Eq. (31), equations (28) and (33) lead 
to 



d x S . Q 
dt T 



(34) 



with Q being the rate of heat exchange between the system and the bath. 
In addition, &\ stands for the entropy production of the system at a given 
value of the external parameter A imposed by the bath. In this case Eq. (27) 
is written as 



-R/k B T 



) 



-AF/k B T 



0(2) 



(35) 



In addition, when the internal irreversible processes in the system have decayed 
(t > t s ), Eq. (34) reduces to the well known Clausius relation SS = 5Q/T , 
valid for quasiestatic processes. 



4 Conclusions 



In this article, we have analyzed the non-equilibrium work relation on the 
grounds of equilibrium thermodynamics and in terms of the underlying micro- 
scopic dynamics based on the BBGKY description. Our analysis is applicable 
no matter how distant from equilibrium the system is. 

We have shown that, in general, in order to be valid non-equilibrium work 
relation's (4) should include a correction taking into account the internal dis- 
sipation in the system. Additionally, we have shown that Eq. (4) corresponds 
to the zero order term in a linear response analysis and thus is merely a near 
equilibrium result (i.e., it applies when the fluctuations are gaussian, in the 
fluctuation-dissipation regime). 
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From the microscopic dynamics of the system, given through the generalized 
Liouville equation and from the entropy postulate (28), we have computed the 
entropy production and provide an estimate of the characteristic time scale 
for the irreversible process taking place in the system, the entropic time t s . 
This quantity depends on the mass, the temperature, the characteristic length 
and the magnitude of the interaction potential and therefore establishes the 
criterion that must be satisfied by a process in order to be quasiestatic. This 
phase-space analysis is the most appropriate for small systems characterized 
by smooth phase-space distribution functions. 

Once again to avoid misunderstandings we emphasize that our microscopic 
analysis is based on a generalization of the Gibbs entropy postulate defined 
through Eq. (28) which is a functional of the distribution vector f(£). This 
distribution vector represents the whole set of s-particle reduced distribution 
functions with s = 0, ,N whose dynamics is given by the generalized Li- 
ouville equation (30). Therefore, our entropy (28) is not a constant of motion 
under the generalized Liouville dynamics (30) as repeatedly has been shown 
in several applications (34; 35; 36; 38). Also incorporated in the theory is the 
presence of a bath which introduces a drive. This has been carried out by 
including the degrees of freedom of the bath into the full phase space through 
the external parameters used to characterize the state of the system. In this 
way, the volumes of the phase space are conserved and the generalized Liou- 
ville equation (33) remains valid. These external parameters originate currents 
entering the generalized Liouville equation and consequently also in the dis- 
tribution function. In this scenario, the non-equilibrium entropy depends on 
the external currents through the distribution functions. 

The present discussion establishes useful quantitative criterion to im- 

plement experiments in good agreement with the theory. 
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Fig. 1. Limiting (solid line) and frequency (dashed line) distributions as a 
functions of the reduced work R* = Rjk^T . The figure shows that (R) ~ AF 
only when the number of experimental measurements N — > oo and the internal 
irreversible processes of the system have relaxed. 
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